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We study ordered states and topological excitations in a quasi- 
two-dimensional magnet, modeled by a square lattice with spins 
s = 1 at all sites, and the Hamiltonian with biquadratic exchange 
interaction between nearest neighbor sites. We propose two 
effective Hamiltonians for description of large-scale excitations 
in the two-dimensional case. They describe excitations of the 
mean field in a nematic phase and a mixed ferromagnetic-nematic 
phase. It is shown that the effective Hamiltonians are minimized 
on configurations with fixed topological charge. These topological 
excitations can arise at low temperatures and cause a destruction 
of a long-range order in the two-dimensional system. 



1. Introduction 

Quasi-two-dimensional magnets have various technolo- 
gical applications. They serve as magnetic films used 
for recording of information, thin ferromagnetic layers 
in Josephson semiconductor junctions, layered resistive 
systems, etc. 

Here, we will not deal with applied aspects of 
theory of magnetism. However, we note that a study 
of two-dimensional systems has also a significant value. 
Studying ordered states, their stability, and excitation 
spectra, we obtain model scenarios of a self-organization 
of a substance with decrease in a temperature or under 
an action of external fields. To support the above- 
presented assertion, it is worth to recall an important 
role played by the Onsager's results [T] on the two- 
dimensional Ising model, or the Kosterlitz-Thouless 
theory of topological phase transitions [213] • This is also 
related to the study of the two-dimensional 0(3)-sigma 
model or a planar Heisenberg magnet [4|. As a natural 
continuation of this trend, we mention a significant 
number of scientific papers devoted to investigations of 
two-dimensional continuous or lattice systems with high 
spins at sites. 

In the present paper we consider a generalized 
Heisenberg magnet taking into account bilinear and 
biquadratic interactions at nearest-neighbor sites of a 



square or cubic lattice. This Hamiltonian for magnets 
of spin s = 1 was proposed and studied long ago 
|6] without any restriction on dimensionality. At the 
beginning of 1970s, an existence of ordered phases 
different from the ferromagnetic or antiferromagnetic 
ones was established by the mean-field methods. In 
particular, if the constant of biquadratic interaction is 
larger than that of bilinear interaction, then a pure 
quadrupole ordering or a spin nematic state can be 
realized in the system [7J[B]. 

It is known that a two-dimensional system with 
a continuous group of symmetries has no long-range 
order at T > (the Mermin-Wagner theorem) . In many 
cases, an instability of ordered phases in two-dimensional 
systems implies an existence of nonlinear topological 
excitations caused by small fluctuations of temperature. 
The role of such excitations in destruction of a long-range 
order is proven within the model of plane rotators 0[3] 
and for the two-dimensional Heisenberg ferromagnet [3| . 

The main result of our work is a proof of an existence 
of topological excitations in the model with biquadratic 
interaction between nearest spins s = 1 at sites of a 
square lattice. We will consider the boundary case of 
a nematic phase, where the constants of bilinear and 
biquadratic interactions are identical. It is known that, 
in this case, energies of both possible phases (nematic 
and ferromagnetic-quadrupole ones) are equal. In order 
to study excitations of nematic phase, we assume that 
K — J — e, and e is a small positive value. It is obvious 
that topological excitations exist at these parameters, 
and differ slightly from those in the case e = 0. If e > 0, 
the manifold of degeneration of a ground state of the 
system is deformed, but the topology is not sensitive 
to smooth deformations. Therefore, our conclusion of 
existence of topological excitations at J — K remains 
valid also in the case K > J > 0. 

The present paper contains two parts. The first part 
is a survey. In the mean field approximation, we reveal 
existence conditions for ordered phases and solve an 
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equation of self-consistent relations for order parameters. 
Comparing with results of other researchers on this 
topic, we obtain conditions of occurrence of a nematic 
state. In the second part, averaging equations of motion 
over coherent states, and passing from a plane square 
lattice to a continuous plane, we obtain formulas for free 
energy of an inhomogeneous distribution of the mean 
field. Topological excitations are described in terms of 
the inhomogeneous distribution. Depending on a choice 
of an equilibrium state and a way of averaging, we 
get two formulas for the free energy: the first formula 
corresponds to excitations of a pure nematic state, and 
the second one is related to excitations in the state 
with nonzero magnetization and quadrupole moment. 
In the case of SU(3)-symmetry, we determine self-dual 
solutions of the problem of minimization. The obtained 
topological excitations give the absolute minimum for 
the free energy, and its value is proportional to a 
topological charge. 

2. The Quantum Model of a Planar Magnet 

Let us consider a plane square lattice containing atoms 
of spin s at each site. Each atom is assigned by 
three spin operators {S^, S^, S 3 } obeying the standard 
commutation relations 

l u n i u n\ — tc ^> n u nm,7 

where the indices a, /3, and 7 run the values {1, 2, 3} 
for each site n, and S nm is the Kronecker symbol. 

Usually, such system is described by the Heisenberg 
Hamiltonian. As s ^ 1 we can include higher orders of the 
exchange interaction in the Hamiltonian. In particular, 
magnets with the biquadratic interaction were studied 
in the 1970s PITO"]. The latter Hamiltonian will be 
considered in what follows. Let 

H = - X){J(S», S n+S ) + K(S n , S n+S ) 2 }, (1) 

where S n denotes the vector (S n , S^, S^) of the spin 
operators at site n, and 5 runs over the nearest-neighbour 
sites. We assume that the exchange integrals J and K 
are positive, i.e. we consider mainly the ferromagnetic 
interaction. 

The operators {S a } are defined over the (2s + 1)- 
dimensional space of an irreducible representation of 
the group SU(2). The spin operators generate the 
complete matrix algebra over this space (the Burnside 
theorem). With respect to the adjoint action adg Q , the 
complete matrix algebra is divided into a direct sum of 



irreducible collections of tensor operators. For example, 
let us consider the case of s = 1. Then for the complete 
matrix algebra over the representation space, we have 
dimMat3 X 3 ~ [9] = [1] + [3] + [5]. Obviously, bases in the 
three- and five-dimensional irreducible collections are 
formed, respectively, by the operators S a , and by the 
tensor operators of weight 2. The latter are the operators 
of quadrupole moment chosen in the form 

Qn 3 = S"S^ + ^ni a 7^ Pi 
Qn' 2] = (Sn) 2 - (Sn) 2 , 

Ql 2 >°J = V3((^) 2 -§). 

A normalization of the operators S a is defined by the 
relation 

(S 1 ) 2 + (S 2 ) 2 + (S 3 ) 2 = s(s + l)I 3 , 

which yields Ti(S a ) 2 = | s (s + l)(2s+ 1). For s = 1 we 
have Tr(S* a ) 2 = 2. We extend such normalization for all 
other operators. 

Now we fix the canonical basis {|+1), |— 1), |0)} in 
the representation space. Then a matrix representation 
of the operators of spin and quadrupole moment is as 
follows: 



1 1 

si = -p o o i , s 2 = — 
^ i i o ^ 



l CP 

& = [ o -l o I , g[ 2 - 0! 

,0 0, 




. Qn 



13 






-A 




i 




( : 


-i 


, Qt 2] = f l 





9 




i / 










Qn 2 



Ql 3 = 

It is easy to see that the proposed matrices are connected 
with the Gell-Mann matrices A a , a = l, 2, . . . , 8, which 
also form a basis in isu(3). The connection is given by 
the linear transformations 

Si = ^(A 4 + A 6 ), S 2 n = ^(A 5 - A 7 ), S n = A 3 , 
Q l2 = k 2l Q[ 2 <°]=A 8 , QW=A U 
Ql? = ^ 2 (h + A7), Q 23 = ^(A 4 -A 6 ). 
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By {P a } we denote the collection of operators {5^, 
Sl Sl QH QH QH g[?' 21 , gL 2 ^ 01 }. It is easy to prove 
that the commutation relations 

hold true. Here, the tensor of structure constants C a b c 
is antisymmetric under a permutation of any pair of 
indices, and its nonzero components are 



123 — ^145 — ^167 — <^264 — <^257 — ^356 
Cl68 = C528 = V3, C437 = 2, 



1. 



In terms of the operators of spin and quadrupole 
moment, Hamiltonian (fTJ) gets the form 



H — — (J — \K) 2_/ « ' "+ 5 



.6 a 



(2) 



,5 a 



where N denotes the total number of sites of the lattice. 
Obviously Hamiltonian remains SU(2)-invariant; hence, 
the operators S 1 " and Q£ are transformed by formulas 
of the adjoint representation 

us%u- 1 = Y,b a ?{u)sl 



UQnU^ 



b 



D ab {U)Ql VU e SU(2), 



where D a P(U) and D ab {U) are matrices of the real 
irreducible representations of SU(2) with dimensions 3 
and 5, respectively. If K = J, then the SU(2)-symmetry 
can be extended to the group SU(3). In this case, the 
Hamiltonian $2$ gets the form 



H 



2 J 2-j z-j n ' 



" na s - \ JN. 



(3) 



.6 a 



To study possible ordered phases of such spin system, 
we use the mean-field approximation. 

3. Mean-Field Approximation. Ordered States 

Now we replace the interaction between spin and 
quadrupole operators that is described by Hamiltonian 
((21 with an effective interaction between the operators 
P£ and the classical mean field. Components of the 
mean field are considered proportional to averages 
(quasiaverages) of the quantum operators {_P°}. The 



Hamiltonian in the mean field approximation has the 
form 



"Hmf = -{J - ^K)^2^2S*(S™ 



+61 



7i, S a 



(4) 



n,5 a 

It is worth to give a warning that a direct calculation 
of the averages and {(Ql)}, for example by 

means of the density matrix p(T)= exp{— w}, results 
in the zero values. This follows from the SU(2)- 
symmetry of Hamiltonian Nonzero values of the 
averages appear if the symmetry is broken. Symmetry 
breaking can be stimulated by an external field, which 
vanishes after specifying an order in the magnetic 
system. The quantities calculated in this way are called 
"quasiaverages" 

Hence, we assume that in our system the nonzero 
quasiaverages and {(Qn)} exist, and form a 

classical 8-component vector field fi a (x n ), a — 1,2, ... ,8. 
In order to obtain nonzero values of {(Ql)}, the external 
field must have nonzero gradient. If the mean field is 
homogeneous, an action of the group SU(2) transforms 
Hamiltonian Q to the form 



HuF^-iJ-^K^SlS 3 )- 



)[2,0] ( Q[2,0] 



)KN 



\KN- 



-J2{(J-lK)S 3 nt i 3 + IKQZ^m) . 

n 

In the case of thermodynamic equilibrium and an 
infinite lattice, the fields /13 = (S 3 ) and /i§ = (Q^ 2 ' ') are 
constant, i.e. have the same values at all points x n 
(a homogeneous mean field). These quantities serve as 
order parameters. Obviously, fj, 3 describes a normalized 
magnetization (a ratio of z-projection of magnetic 
moment to a saturation magnetization), and /i§ is 
analogously related to a quadrupole moment. 

In the mean field approximation, it is easy to 
calculate a partition function for the homogeneous case 



NZ([i 3 ,[i 8 ,T) = Tre 



Tr e 



where /imf denotes a one-site Hamiltonian 
h MF = -(J - kK)^S 3 - \K^Q™ 



\K. (5) 



2 2 '-"o^ 3 

The introduced mean field makes sense if self- 
consistent relations are held: 

M3 = \S )MF 



Tr e 
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Ms = (Q [2 ' 01 )a/f = 



TrQ[ 2 >°] 



Tr e 



Nh MF 



these relations serve as an analog of the Weiss 
equation from theory of ferromagnetism. The averages 
of operators can be presented via the partition function: 



M3 



kT d 
(J-f)5^ 



lnZ(iJ, 3 ,fj, 8 ,T), 



2kT 9 1 7( T\ 
M8 = 77 o hlZ(^3,^ 8 ,T). 

For the system described by one-site Hamiltonian ([3]), 
the self-consistent relations get the form 



2sh 



M3 



Ms 



(J-f )M3 



exp 



\ 2fcT J 



2ch 



fcT 



ch 



(J-f )w 

kT 



exp 



\ 2fcT J 



^ex P {-^} 



2ch 



(6) 



fcT 



Note, that the true averages are always less than their 
expectation values calculated from the self-consistent 
relations. Therefore, solutions of © have a qualitative 
sense only. 

Here we analyze Eqs. ((HI and make comparison 
with results described in the literature. The obvious 
solution corresponds to the paramagnetic state with 
jL«3 = and fig — 0; this state is realized at temperatures 
kT > J — K/2. At the same time, this inequality shows 
that the model, oriented to ferromagnetic materials, 
gives a meaningful result only in the region J — K/2 < 0, 
which contains areas with the ferromagnetic and 
quadrupole orderings, according to the well-known phase 
diagram (Fig. 1) for the bilinear-biquadratic s = l 1- 
dimensional spin model [12]. 




Fig. 1. Phase diagram of a one-dimensional system of spins s = 1 

In the case of K < 0, the self-consistent relations 
have a unique nontrivial solution, corresponding to the 
ferromagnetic ordering, because this solution tends to 
the boundary values /i3 = 1 and A*8 = ^ as temperature 
decreases to zero. The critical temperature of transition 
from a ferromagnetic state into a paramagnetic one 
is determined in terms of the constants J and K as 
T c = &(J-K/2). 

In the case of K > (the light-grey region in 
Fig. 1), Eqs. ([6]) have more than one nontrivial solution: 
two solutions corresponding to ferromagnetic states 
with the boundary values = 1 and /i 3 2 ^ =2/3 (and 
the corresponding values of /ig), and two solutions 
describing nematic states (/13 = 0) with the boundary 
values /Xg = ^| and /i^ — Existence of four ordered 
states in ferromagnets is also reported in [7]: they 
are a ferromagnetic state with /ig = s, a quadrupole 
(or nematic) state with /i3 = 0, /ig 1 '' = — s(s + 1)/V3, a 
partially ordered quadrupole state with /i 3 = 0, > 0, 

and a partially ordered ferromagnetic state with fj,^ < s. 
Partially ordered states are unstable [7]. 

Analyzing the temperature evolution of solutions of 
© as K > 0, J > 0, we reveal two critical temperatures, 
which are solutions of the equation 



J - K/2 
kT 



- l) = exp I 



K / 
kf V ~ 2(7 



3kT 



K/2) 



An obvious solution is T cl = Jr(J — K/2). The other 
solution T C 2 is calculated numerically. In the region 
J> K, i.e. for ferromagnetic materials, the temperature 
T C 2 is less than T c i, whereas the reversed situation 
takes place for nematics in the region J < K. At smaller 
critical temperature the solution disappears. Then 
only the solution fi^ exists in a certain interval of 
temperatures. A comparison with results of the paper [8] 
shows that at a higher critical temperature we have a 
second order phase transition from a ferromagnetic state 
into a paramagnetic one. 

A somewhat different behavior is observed 
in materials with J w K, or more precisely 
0.67 < J/K < 1.22, i.e. on the boundary between the 
ferromagnetic and nematic regions (the dark-grey 

region in Fig. 1). Disappearing at a lower critical 

(2) 

temperature, the solution /i 3 arises again at a higher 
critical temperature, and grows continuously from zero 
toward . When two solutions coincide at a certain 
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temperature To, they disappear by jump. We may 
conclude that the phase transition from a ferromagnetic 
state to a paramagnetic one is a transition of the first 
order. This well agrees with results of the paper |5], 
where the change of a second order phase transition 
into a first order one in the region 2 J/3 < K < J was 
considered (for ferromagnetic materials). 

Further we consider the case J = K, corresponding to 
the boundary between the ferromagnetic and quadrupole 
regions. Moreover, as J = K the quantum Hamiltonian 

and the Hamiltonian in the mean-field approximation 
are SU(3)-invariant. The latter has the form 



H 



MF 



(7) 



4. Equations of Motion for Large-Scale 
Fluctuations of the Mean Field 

Now we return to the quantum SU(3)-invariant spin 
model with Hamiltonian Q. The Heisenberg evolution 
equations for the operators P° have the form 



dP a - 



(8) 



We assuming that the system is ordered state, and 
take an average of the both sides of Eq. ([5]) over 
Heisenberg (time independent) coherent states 

W»)> - -f=(ci(n)|l> + c_i(n)|-l) + co(n)|0», 



of the discrete lattice. Such transition is well known 
for an SU(2)-magnet [2] and underlies the macroscopic 
phenomenological theory of magnetism [15]. In the 
continuous theory dynamical variables are densities of 
averaged spin and quadrupole moments: 

M a (x) = lim — Ha(n)5 X:Xn = V] fjt a (n)S(x-x n ). 



nes 



neS 



Here, S is a 'physically' infinitesimal region of the lattice, 
and 5(x—x n ) is the Dirac delta-function, which has the 
dimension of reciprocal area. A Poisson bracket on the 
space of {M a (x)} is calculated by the formula 

{M a (x),M b (y)} = C abc M c (x)5(x - y). 

In what follows, we deal with the dimensionless field 
/i a (x) = I 2 M a (x) , where I is a distance between the 
nearest neighbors of the square lattice. 

Considering (j±l,fc), (j,k±l) as the nearest 
neighbors of the site n — (j, k), in Eqs. (fTU| we perform a 
transition from the discrete variable x n to a continuous 
one x. Then the field {fi a (x)} satisfies the equations 

,dn a (x) 5H cS 
ft = ?4ff} = -CabcUbs , (11) 



dt 



SfJ-c 



where 



A suitable representation of the system of Hamilton 
equations pip is the matrix equation 



\ Cl \ 2 + \c^\ 2 + \c \ 2 = l. 



On the other hand, an averaging can be performed 
by means of the density matrix (thermodynamical 
averaging) as T > [13]. In both cases, we suppose 



(CO « {Pn)(PD = Ma(n)/x 6 (m), 



(9) 



dfi 
Ik 



2J1 2 



(12) 



where 



f 1 = 



2 ^ 



[ i a (x)P a . 



i.e. we neglect correlations between fluctuations of the 
quantum fields P£ . Then we obtain the following system 
of Hamiltonian equations for the averages /i a (n): 

n — - — = C abc n b (n)- — — = {n a {n), {%}}. (10) 



dt 



dfic(n) 



Taking © into account we have (H) = (Hmf)- 

In order to investigate large-scale fluctuations of the 
field fi a (n) we take a two-dimensional continuum instead 



Obviously, /t is a Hermitian 3x3 matrix, the bracket [•, •] 
denotes a matrix commutator, and A = + g?2 is the 
Laplace operator over a two-dimensional space. Equation 
(fTTj) generalizes the Landau-Lifshits equation [TB] for 
isotropic ferromagnets to the case of 8-component mean 
field /J, a (x). The Landau-Lifshits equation is well-known 
in the macroscopic theory of magnetism, and suitable for 
describing large-scale excitations in planar magnets and 
exploring the ferromagnetic resonance. 
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5. Free Energy and Topological Charge 

As mentioned above, a Hamiltonian coincides with a free 
energy at constant temperature and volume. Therefore, 
we use the notion of free energy in what follows. As 
shown in Appendix (Section 6) by means of an algebraic 
approach, Eq. (j!2[) coincides with the two-dimensional 
generalization of Eq. (f2"5"j) on a degenerate orbit. This 
equation is a Hamiltonian one, though it is nonintegrable 
in the two-dimensional case; and its Hamiltonian can be 
used as an effective free energy for the spin system in 
question, namely: 



•eff 



2 

3h 



A,y) 



dxdy. 



Obviously, Ff s is a part of the total free energy of a 
magnet, and arises from an inhomogeneous distribution 
of the average values {/i a (x)}. 

The algebraic approach yields one more equation of 
motion, associated with a generic orbit. Evidently, this 
equation can be obtained from the quantum-mechanical 



approach, by performing a relevant averaging that takes 
correlations into account. Therefore, we consider another 
effective free energy 



1 



2(3/ 2 



hi) 



(^Cj(mL 



Ma,, 



T a>v ) + 3h (Tl x + Tl y ))dxdy. 

Equations of extremals for the functionals of free energy 
are the two-dimensional generalization of Eqs. (1251) and 

If equations of constrains determining orbits are 
solved or, in other words, orbits are parameterized, then 
the formulas for the free energy can be simplified. Orbits 
of co-adjoint representation of semisimple compact Lie 
groups are compact Kahlerian manifolds. Therefore, 
it is suitable to use a complex parameterization. In 
order to parameterize a generic orbit, it requires three 
complex variables ui\, w 2 , d 11)3. Explicit formulas for the 
parameterization of a generic orbit are the following: 



Mi = 



M2 = 



M3 



M4 



Me 



M7 



Ms 



- y/3q W2 + w 3 + w 2 + u>3 m (1 — wi)(w 3 — W1W2) + (1 — wi)(w 3 — wiw 2 ) 

y/2 1 + |wi| 2 + \w 3 - wiw 2 \ 2 

im (1 + u>i)(w 3 — W1W2) — (1 + wi)(w 3 — W1W2) 



2V2 

m — V3q 



1 + \w 2 \ 2 + \w 3 \ 2 

W3 - W 2 - W3 + W 2 



2iy/2 1 
m — V3q 



\w 2 



2 



m 



1 + |l0 2 | 2 + |t0 3 | 
W2W3 - W2W3 



2i 

to — V3q 
2^2 



12 



W3 — W2 + W3 — W 2 



y/2 1 + \w\\ 2 + \w 3 - WiW 2 | 2 

m(l-| Wl | 2 ) 

1 + |wi| 2 + \W3 - WiW 2 | 2 ' 

im{w\ — w\) 

+ \wi\ 2 + \W3 — W1W2I 2 ' 

m (1 + toi ) (703 - W1W2) + (1 + wi)(w 3 - i0ii0 2 ) 



1 



to — v3 q W2 
2iV2 1 
to — \/3q 



\w 2 \ 
W3-W2 

12 



r^3 



V2 



1 + |l0i| 2 + \w 3 - t0i W 2 | 2 



W2\ + \ w 3 
W2W3 + W 2 W3 



W 3 
2 



im (1 - Wi)(w 3 ~ Wiw 2 ) - (1 - Wi)(w 3 - w 1 w 2 ) 
\pl 1 + |wi| 2 + \w 3 - wiw 2 \ 2 

m(wi + w\) 



1 



102 



™3 



\/3< 



w 2 



2V3 



1 



\w 2 \ 



TO 

71 



103 - WiW 2 \ 

2 -2\w3 



101 



101102 



1 + |l0l| 2 + \w 3 - I0il0 2 I 2 



(13) 



Here, m and q denote boundary values of the quantities the expressions for free energy get the form 
[13 and fig, respectively. For a degenerate orbit, it is 

sufficient to have two complex variables, for instance w 2 -pes _ f - [ ^ Wa _|_ ^ Wa j (14) 

and W3; in this case, we assign w± = 0. J „ a ^ \ dz dz dz dz J 



After the restriction onto an orbit by formulas (|13[) . where g a a denote components of a metrics on the orbit 
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and real coordinates x, y on a plane are changed into 
complex ones z, z. 

While coadjoint orbits are Kahlerial manifolds they 
possess Kahler potentials, which generate related metrics 
tensor g and 2-form h; their components are calculated 
by the formulas 



<9 2 $ 



<9 2 $ 



y ° p dw a dwp' Q/3 dw a dw p ' 
A 2-form gives rise to a topological charge 

1 f ^ fdw a dwp dw a dwp\ 

J a/3 v 7 

which means a degree of mapping of a plane into an 
orbit, realized by the function w\, W2, W3. 

On a degenerate orbit of SU(3) the function 
$2 = hi(l + |w2| 2 + |w3| 2 ) serves as a Kahlerian 
potential, and the metrics tensor from (fT¥|) is a Kahlerian 
one. Then the formulas for topological charge and free 
energy differ only in a sign ('+' for a free energy, and 
' — ' for a topological charge), hence, 

m >4tt|Q|. 

The equality holds if the second term in the brackets 
vanishes, that happens if the functions {w a } are 
holomorphic or antiholomorphic. 

Consequently, holomorphic functions form a class 
of solutions with quadrupole ordering (to = 0) that 
correspond to the minimum of free energy. The same 
takes place for antiholomorphic functions. 

Now we consider the case of a generic orbit. 
The cohomology class of rank 2 for the orbit 
is two-dimensional, then there exist two basis 2- 
forms, generated by the Kahlerian potentials $ 2 > 
and $1 = ln(l + |wi| 2 + \ws — Wiu^l 2 ). As a unique 
Kahlerian potential we take the one corresponding to 
the Kirillov-Costant-Suriau form 



<f> = TO$i $ 2 - 



Then all conclusions relative to a degenerate orbit 
can be extended to a generic one. That is, the class of 
solutions with ferromagnetic ordering that correspond to 
the minimum of free energy are realized by holomorphic 
(or antiholomorphic) functions. 

5.1. Large-scale topological excitations 

Excitations of a state with quadrupole ordering are 

described by spatially inhomogeneous distributions 

of the 8-component vector field /x a (a;), living on a 
degenerate orbit 



O(fX 3 = 0,» 8 ) 



SU(3) 



SU(2) x U(l) 



Mappings of topological charge 2 can be modeled by the 
holomorphic functions 



w 2 (z) 



Z\ 



W3(z) 



a 2 



Z2 



where 01, a 2 , Z\, and z 2 are fixed complex numbers. 

The components /i 3 and fj,g of the mean field, 
whose boundary values are called order parameters, are 
presented in Figs. 2 and 3 (we show a cut along the 
straight line joining poles of the functions w%(z) and 
W3(z), the origin of coordinates being at a middle of 
the interval ziz 2 ). In Figs. 2 and 3, q is a value of the 
component fig at an initial point of a degenerate orbit 
(M3 = 0). 




(15) 



Fig. 2. Contour of fi:j(z, 2), ^3(00) — > 



Generally speaking, the metrics tensor of free energy (Tnj 
is not Kahlerian. However, we can construct a Kahlerian 
metrics of the form 



x-eff _ 

•'2 — 



1 



T a<y ) + C 3 (Tl x + Tl y ))dxdy, 
by choosing appropriate coefficients C±, C 2 , C3. 



q 




\ 2 


1 / 


\ Zi -\ z 




/ q 

2 
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Fig. 3. Contour of fJ,s(z. z), ^tg(oo) — > q 

These excitations are analogous to Belavin-Polyakov 
solitons, well-known for planar Heisenberg ferromagnets 
(the quantities a\ and a 2 define widths of solitons, 
and the quantities Z\, z 2 give their positions). It is 
easy to see that an energy of a configuration does 
not depend on a width of soliton, that proves scale 
invariance of the energy in the two-dimensional case. 
Hence, topological perturbations can have arbitrary 
large size. Such instability (an unrestricted increase of 
soliton width without pumping of energy) can cause a 
destruction of the nematic order in the considered model. 



6. Appendix. Integrability of SU(3)-Symmetric 
Equations of the Landau Lifshits Type in a 
One-Dimensional Space 

It is known that system of equations (TT3]) in the one- 
dimensional case is an integrable Hamiltonian system on 
a degenerate orbit of the group SU(3) [T7]. We generalize 
Eq. (fT2|) to the case of a generic orbit. The algebraic- 
geometric nature of equations like (fT2j) is revealed in 
the frame of so called orbit approach to nonlinear 
Hamiltonian systems. Below, we construct integrable 
Hamiltonian equations on orbits of the group SU(3). 

Consider polynomials in a complex variable A, 
whose coefficients are anti-Hermitian matrices of the 
algebra su(3). We denote the set of polynomials by 
g + ~ su(3) ® V{\), where V{\) is a ring of polynomials 
with the standard multiplication operation. If A. B G g + 
have the form A(X) = J2 n ^»A", B(X) = £ fc B k X k , then 



[A, 5]=]T[i„, B k ]\ n+k eB+- 



(16) 



■ii.k 



Operation (fTBj) shows a structure of graded Lie algebra in 
g + . Let X" = X n X a be a basis in g + , where X a = — |A , 
a = 1, 2, . . . , 8, {A a }® =1 denote the Gell-Mann matrices. 
In we introduce a bilinear ad-invariant form 



(A, B) = — 2resA TrA(A)i?(A), 



(17) 



where N + 1 is the maximum degree of matrix 
polynomials A and B. Let M. — (5+)* be a space dual to 
g + with respect to form (fTTj) . A collection of the linear 
forms 



JV 



^ = E E + (&*3 + A 



N+l 



creates a closed ad-invariant subspace M N in M.. The 
coordinates ^™ of ^(A)^^^ are calculated by the 
formula 

eH£(A),A-" +Ar+i ). 

In the linear space M. N with coordinates 
n = 0, 1,...,N, we define the Lie-Poisson bracket 



{a,/ 2 } = EE^ 



, dfi df 2 

ra,n a.b u 



(18) 



with the Poisson tensor field 



WT = <£(A), [x~ m+N+1 , x b - n+N+1 ]). 



We also define two ad-invariant functions /2(A) and /3(A) 
by the formulas 

/ 2 (A) = -2Tr£ 2 (A) = ^(A), 

a 

h(X) = -4zTrf(A) = «U:€«(A)6(A)€c(A), 

Here, d abc = -2i Tr(X a X b X c + X b X a X c ), and $„(A) is a 
polynomial 

ea(A) = e Q ° + ^A + ^A 2 + --- + e +1 A Ar+1 . 

The invariant functions are also polynomials in the 
complex parameter A: 

7 2 (A) = / l0 + ^ 1 A + --- + / l2Ar+2 A 2A ' +2 , 

/ 3 (A) = jo + /lA + • • • + f;x.:;\ :iX \ 

It is easy to prove that the coefficients ho, ■ ■ ■ , /ijv+i, 
/o, ... Jn+i are annihilators of bracket (I18p . Fixing 
them we obtain the system of algebraic equations 

h n = const, /„ = const, n = 0, . . . , N + 1, (19) 

which give an embedding of the orbit O n+1 of dimension 
6(7V + 1) into the linear space M N+1 . The rest of 
coefficients {h N+2 , . . . , h 2 N+2, In +2, • • ■ , /3W+3} form 
a pairwise commutative collection of integrals of motion, 
which is necessary to integrate the Hamiltonian system. 
We are interested in the functions /ia/+ 2 and /ijv+3 and 
the corresponding Hamiltonian equations. In particular, 
the Hamiltonian tiN +2 gives rise to so-called stationary 
equations. In terms of the coordinates £™, they are 



dx 



9 f tOen+l 

£.1 abcC.c^b ' 



(20) 



n— a— 1 



where {f a bc} are antisymmetric structure constants of 
the algebra su(3) in the basis of Gell-Mann matrices: 

[X a , Xb] — fabcX c , 
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/l23 — 1) /458 — /786 — ^} , 

/l47 = /l65 = ^246 = /257 = /345 = /376 = g> 

By a; we denote a "time" with respect to the Hamiltonian 
hN+2, which corresponds to evolution equations 



dt 



(21) 



Equations (12"01) and ([2"Tj) are consistent, 
because the corresponding Hamiltonians commute: 
{h-N+2, ^7V+3} — 0. Hence, we can determine evolution 
(|2"Tj) on trajectories of system (|20| . i.e. we suppose the 
dynamical variables £™ in Eq. (j2"Tj) to be dependent on 
x. Combining ([2U)l and (|2"Tj) , we have 



da 



dt - 2 f<*°£$- >)r 



(22) 



The variables can be expressed in terms of 

and then (f2"2"|) becomes a closed system of partial 

differential equations for For it is necessary 

to solve the degenerate system of equations 



^ — ^Jabc^c^b ' 



(23) 



that becomes possible after restriction to an orbit 
6.1. Classification of orbits 

It follows from Eqs. that the orbit O n+1 in M N+1 
has the structure of a vector bundle over a co-adjoint 
orbit of the group SU(3). Hence, a classification of the 
orbits O n+1 is reduced to that of orbits of SU(3). 

Since the group SU(3) is simple, we have g*=g. 
Therefore, we consider also as coordinates in the 
algebra g ~ su(3). Then a generic element £ £ su(3) is 
represented by the matrix 




-& + -71^8 

e 6 ° + ^7° 



e 6 °-^7° 
£ 8 



(24) 



_2 

x/3 



Let £(0) be a fixed element of su(3). By definition, 
the set O ? - (0) ={g|(0)5 _1 , VgeSU(3)} is an orbit of 

SU(3) through the point £(0). Elements g' of the group 
SU(3) such that g'i^g'^ 1 =|(0) form the stationary 
subgroup at the point £(0). An orbit O^, y being a 
homogeneous space, is a coset space SU(3)/G'~O^ s, 
where G' denotes the stationary subgroup. 



A maximal commutative subalgebra of a semisimple 
algebra g, which is called Cartan, can be diagonalized. 
The Cartan subalgebra of su(3) is formed from diagonal 
matrices depending on two parameters £3 and £g. It is 
well known that a proper transformation puts any 
element £ £ g into the Cartan subalgebra of g. This yields 
that each orbit intersects the Cartan subalgebra at least 
once. In fact, there is more than one intersection point 
number, more precisely as many as an order of the Weyl 
group W(G). We discuss this in what follows. 

Nontrivial similarity transformations £— > gig^ 1 that 
preserve a subalgebra f) form a discrete subgroup 
W(G)cG, which is called a Weyl group [TH]. An 
action of the group W(G) on the subalgebra f) = f)* 
is generated by reflections in planes orthogonal to 
simple roots. The Weyl group of SU(3) is generated 
by two reflections o~\ and (72 in the planes shown in 
Fig. 4 by dotted lines. The full Weyl group consists of 
six elements {e, o~\, 02, uia^, <72ai, o\O20~\ — 0-2010-2}, 
and is isomorphic to the group of permutations S3, 
ordS 3 =3!. 

Weyl chamber 



a 2 








'' a i\/ 






0-, 



Fig. 4. Root diagram of the group SU(3) 

Since a Weyl group acts over a Cartan subalgebra, 
every point cr^(0)cr _1 , a €W(G), belongs to an orbit 
through £(0) € fj. The group W(G) acts efficiently 
(change an element £ into another, if the element does 
not belong to reflection planes). An open domain in a 
Cartan subalgebra where a Weyl group acts efficiently, is 
called a Weyl chamber (see Fig. 4) . Elements of different 
Weyl chambers are adjoint by elements aGW(G), this 
implies that each orbit through a point £(0) of Weyl 
chamber intersects a Cartan subalgebra as many times 
as an order of W(G). If £(0) is an interior point point 
of a Weyl chamber, we call the orbit a generic one, and 
call the points cr£(0)cr -1 , Ver G W(G) poles of an orbit. 

For a generic orbit a stationary subgroup G' 
coincides with a maximum torus T r for a group G 
of rank r. In the case of group SU(3), we have 
T 2 =U(1) x U(l). Hence, generic orbits are coset spaces 
O gen ~SU(3)/U(l)xU(l). 
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If an initial point £(0) belongs to a wall of a Weyl 
chamber (in the case of SU(3), belongs to one of 
the reflection lines), then we deal with a degenerate 
orbit. In this case, a stationary subgroup G' contains 
a semisimple subgroup generated by roots orthogonal 
to an initial point £(0). Consider the group SU(3). 
If £(0) lies the vertical reflection line, aj and —a.\ 
are orthogonal to this element. The corresponding 
sl(2)-triple {X ai , X_ ai , H ai = [X ai , X_ ai ]} generates 
a subgroup SU(2) C SU(3). Obviously, the element 
£(0) = ^= £g diag(l, 1, — 2) is invariant under a 

transformation g'£(0)g'~ , where g' is the unitary matrix 



a 


P 






a* 


( 


\ o 








J<p/2 






Jtp/2 








—itp 



Hence, g 1 G SU(2) x U(l), and a degenerate orbit is a 
coset space O dog ~ SU(3)/SU(2) x U(l). 

6.2. Equations on orbits and their Hamiltonians 

Return to construction of integrable systems on orbits 
of loop groups, and consider Eq. ((2"2"j) . In order to solve 
this equation we have to restrict the degenerate system 
(|2"3"|) into an orbit and solve. 

If £° G Cdeg, then the matrix 2f abc Q has rank 4, and 
its inversion gives a solution 

pi _ _2_ f cO^O j hi cO 
— 3h JabcSb'scx ~ 2h Sai 

where the constants ho, hi define an orbit by Eqs. ([TO]) . 

If £° £ Cgon, then the matrix 2f abc Q has rank 6 and 
is invertable on a generic orbit. Then we have 

+3h f abc T]%rf cx - 6fofabc£bVix) + 

, 2/ /i-3hgfti >0 , 3/0^1-2/iq/! 
+ 6(/^-hg) «« + 6(/^-^) 

where 77^ = rfabc£b£c > an( i the constants /iq> /o> /1 
come from Eqs. (fT9|) . 

Substituting the obtained expressions in the right- 
hand side of (f2"^j) , we get two equations for the functions 
Mx,t) = £: 



dt 3h ° 

dtU 

dt 2< - h l 



3_3j2^ {h(jf abc£,b£,c,xx ^fofabcS,bVc,x 



-3A 



ofabcVbVc^x — 3fofabc11b£c,xxj~ 



I 2/ /i-3fe fei £ , 3/ fei-2fe /i Pc.fi 



Let /ii = in Eq. (f?5|) , and replace the variables £ a by 
/i a . Then its generalization to the two-dimensional case 
gets the form 



dt 



where A = + gp- • Obviously, this equation has the 
Hamiltonian 



u 



eff 



1 

12/10 



«) dxdy. 



It is easy to see that (1231) coincides with a one- 
dimensional analog of (fT!?)) . In other words, Eq. (jT2"|) can 
be considered as a two-dimensional generalization of the 
integrable equation (l23|l . 

In the same way we treat with Eq. ()26[) . namely, 
replace £ Q by fi a and assign fi =h\ =0. Generalized to 
two dimensions, the obtained equations get the form 



dt 



8(fcg-3/ ( 



F2) (h-oCabcfJ'bA-lJ.c — if Q C abc rj b A^ c 



+3h C abc r] b Ari c - 3foC abc n b Ar] c 



(27) 



Here, ?7 a are quadratic forms in fi a : rj a — d abc [i b \x ci 
where d abc = \ Tr(P a P b P c + P b P a P c ). Obviously, Eq. 
(|2T|) is Hamiltonian, and give the following effective 
Hamiltonian 



n 



eff 



16(/ig-3/ 2 ) jfi^O^x + h lA,,y~ 

-6foVa,xVa,x ~ 6foHa,yVa, y + 3h r% x + 3h n r^ y ) dxdy. 
7. Conclusions 

In the present paper, we have constructed nonlinear 
stationary excitations appearing in the nematic phase 
of a planar magnet of spin s = 1, modeled by a square 
lattice with a biquadratic interaction between nearest- 
neighbor sites. These excitations are characterized by 
an integer topological charge, and reveal themselves 
as regions with nonzero magnetization and mean 
quadrupole moment. Topological excitations in a two- 
dimensional system (without taking into account an 
anisotropy and a demagnetizing field) can increase 
unrestrictedly without pumping of energy. This destroys 
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a nematic state in the system, according to the Mermin- 
Wagner theorem on absence of a long-range order in one- 
and two-dimensional systems. 
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BnOP5LHKOBAHI CTAHH TA HEJIIHIHHI 
BEJIHKOMACIHTABHI 3BY^>KEHHH 

y njiocKOMy mafhethky 

31 CniHOM s = l 

K).M. BepnaupKa, II. I. ToAod 

P e 3 K) m e 

7tocjiifl>KeHO BnopsiflKOBam CTaHii Ta Tonojiorinm 36yjx>KeHHH 
y KBa3iflBOBHMipHOMy MaraeTHKy, 3MO^ejibOBaHOMy KBajxpaTHOio 
f paTKOio 3i cnmaMii s = 1 y By3jiax Ta raMijiBTOHiaHOM 3 SiKBafl- 
paTHOio o6MiHHOK> B3aeMOflieio Hait6jin>KHHx cycifliB. 3anponoHO- 
Baeo jxsa ecbeKTHBHiix raMijiBTOHiaHH jjjih oniicy BejiHKOMacniTaS- 
hhx 36ys*ceHb y CTporo flBOBHMipHOMy BHnaflKy. OflHH 3 hhx onn- 
cye 36yfljKeHHH cepeflHBoro nojia b hcmcitemhih <pa3i, minim — y 
3MiiuamH cbepoMarHiTHO-HeMaTiiHHiii cba3i. noKa3aHO, m,o edpeK- 
THBHi raMijiBTOHiaHH MiHiivii3yiOTbC5i Ha KOHcbirypaii,i5ix, Hid MaiOTB 
cpiKCOBaHHH TonojioriHHHii 3apHfl. IJ,i TonojioriiHi 36yfl>KeHHH mo- 
»yTB BHHHKaTH npH He3HaiHHx TeivmepaTypax i 6yTH npHHHHOio 
pyHHyBaHHH flajiBHBoro nopHflKy y CTporo flBOBHMipHiii cncTeMi. 
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